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Abstract 

The asymmetric shape of the nebula around //-Carinae (Homunculus) can be ex- 
plained by a spherical expansion in a non-homogeneous medium. Two models are 
analyzed: an exponential and an inverse power law dependence for the density as a 
function of distance from the equatorial plane. The presence of a medium with vari- 
able density along the polar direction progressively converts the original spherical 
shell into a bipolar nebula. In the case of the nebula around r/-Carinae, we know the 
time elapsed since the great outburst in 1840. An exact match between observed 
radii and velocities can be obtained by fine tuning the parameters involved, such as 
initial radius, initial velocity and the typical scale that characterizes the gradient in 
density. The observed radius and velocity of the Homunculus as a function of the 
polar angle in spherical coordinates can be compared with the corresponding simu- 
lated data by introducing the efficiency in a single or multiple directions. Once the 
3D spatial structure of the Homunculus is obtained, we can compose the image by 
integrating along the line of sight. In order to simulate the observed image, we have 
considered a bipolar nebula with constant thickness and an optically thin emitting 
layer. Some simulated cuts of the relative intensity are reported and may represent 
a useful reference for the astronomical cuts. 

Key words: ISM: Molecules, Stars: individual: eta- Carinae, Stars: Mass Loss, 
Stars: Winds, Outflows 



1 Introduction 



The nebula around rj- Carinae was discovered by iThackerayl (119491 ) and the 
name "the Homunculus" arises from the fact that on the photographic plates 
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it resembled a small plump man, see iGaviolal (119501) . M ore details on the 



various aspects of 77-Carinae can be found in ISmithl (120091 ). The structure of 
the Homunculus Nebula around 77-Carinae has been analyzed with different 
models, we cite some of them: 



The shape and kinematics is explained by the interaction of the winds ex- 
pelled by the central star at different injection velocities, see lickl fll988f ). 
The possibility that the nebulae around luminous blu e variables ( LBVs ) 
are shaped by interacting winds has been analyzed by Not a et aD ( 1995 ). 
In this case a density contrast profile of the form p = p (l + 5 cos 4 6) was 
used where is the angle to the equatorial plane. 

The origin and evolution of the bipolar nebula has been modeled by a nu- 
merical two-dimensional gasdynamic model wher e a stellar w i nd int eracts 
with an aspherical circumstellar environment, see iFrank et al.l (119951 ). 
Cooling models form ballistic flows (that is, a pair of cones each with a 
sph erical base) whose la t eral e dges become wrinkled by shear instabilities, 
sec 



Dwarkadas fc Balickl (119981 ). 



The scaling relations derived from the theory of radiatively driven winds 
can model the outflows from luminous blue variable (LBV) stars, taking 
account of stellar rotation and the associated latitudinal variation of the 
stellar flux due to gravity darkening. In particular for a star rotating close 
to its critical speed, the decrease in effective gravity near the equator and 
the associated decrease in the equatorial wind speed results naturally in a 
bipolar, prolate intera c tion f ront, and therefore in an asymmetric wind, see 



Dwarkadas k Owockil ((2002). 



Two oppositely ejected jets inflate two lobes (or b ubbles ) repr e sentin g a 



Sokerl (12004120071 ). 



unified model for the formation of bipolar lobes, see 
A two-dime nsional, time-de p enden t hydrodynamical simulation of radiative 
cooling, see iGonzalez et al.l (120041 ) . 

Launch of material normal to the surface of the oblate rotating star with an 
init ial kick velocity tha t scale s approximately with the local escape speed, 



sec 



Smith fc Townsendl (I2007f ). 



A 3D m odel of wind-wind co llision for X-ray emission from a supermassive 



star, see 



Parkin et all (120091 ). 



Two-dimensional hydrodynamical simulations of the eruptiv e events of the 



1840s ( the great outburst) and 1890s (the minor outburst), see lGonzalez et al. 
(|2010ah . 



The models cited leave some questions unanswered or only partially answered: 

• Which is the law of motion which regulates the expansion? 

• Is it possible to model the complex three-dimensional (3D) behaviour of the 
velocity field of the expanding nebula? 

• Is it possible to make an evaluation of the reliability of the numerical results 
on radius and velocity compared to observed values? 
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• Is it possible to evaluate the intensity of the H 2 image of the nebula? 

• Is it possible to build cuts of the model intensity which can be compared 
with existing observations? 

In order to answer these questions, Section 2 describes three observed mor- 
phologies of PNs, Section 3 analyzes two different laws of motion with their 
associated velocities which model the aspherical expansion and Section 4 con- 
tains detailed information on how to build the image of the Homunculus as 
well as some sectional cuts through the relative intensity of emission. 



2 Homunculus properties 



This Section reviews the basic data of the Homunculus as well as an ad 
hoc formula for the latitude-dependent wind. The star ^-Carinae had a great 
outburst in 1840 and at the moment of writing p resents a bipo lar shape 
called the Homunculus, its distance is 2250 pc, see ISmith I ( 120021 ). A more 
refi ned classificatio n disti nguishes between the large and little Homunculus, 



sec 



Gonzalez et al.l ( 120061 ). The Homunculus has been observed at d i fferent 



wavele ngths suc h as the ultravio l et and infrared bv lSmith et al. (2004); ISmith! 
(l2009h . x -ray by[Corcoran et al.l (l2004ah . [Fei7l A16435 bvlSmithl (12005^ am- 
monia bv lSmith et al. I (l2006h , ra dio-continuum by lGonzalez et al.l (I2006J) , near- 
infrared by lTeodoro et al.l (120081 ) and scandium and chromium lines by lBautista et al 
(120091 ). 



Referring to Table 1 in ISmith! (120061 ). we can fix the major radius at 



22014 AU (0.106 pc) and the equatorial radius at 2100 AU (0.01 pc). The 
expansion velocity rises from rj 93 km/s at th e equator to fa 648 km/s in 
the polar direction, see Table 1 and Figure 4 in ISmith! (120061). The thi ckness 
of the H 2 shell is roughly 2 — 3% of the polar radius, see ISmithl (120061 ). The 
expansion speed of the outer H 2 shell of ?7-Carinae has been fitted with the 
following latitude dependent velocity as given by 



,2 A cos(2 6) 



Vl (J +e 2Acos(2G)) 



where the parameter A controls the shape of the Homunculus, G is the polar 
angle; v-\ and Vo. are the velocities in the polar and equatorial d irection, see 

1 reports the data of ISmith! ( 120061 ) as well the 



Gonzalez et al. ( 


2010a 


). Figure 


fit of 


Gonzalez et al. ( 


2010a). 
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Fig. 1. Expansion velocity versus latitude of the H2 shell (Smith 2006, open stars) 
and fit as given by for mula (1) (full line). The fitting parameters are the same as 
Gonzalez et al.l (|2010al ): A=1.9, v\ = 670 km/s and V2 = 100 km/s. 



3 Law of motion 



This Section presents a standard solution for the law of motion in an inter- 
stellar medium with constant density. The assumption of constant density can 
be generalized considering a medium which presents a decreasing behaviour 
in density with increasing distance from the main plane of symmetry: in par- 
ticular we considered an exponential and a power law decrease in density. As 
a consequence two new laws of motion are deduced. In all three cases, the 
velocity is deduced assuming the conservation of momentum and has a simple 
analytical form. The resulting first order equations for the radius have sep- 
arable variables and therefore the solution can be found. In the first case is 
possible to derive an analytical expression for the radius and in the two other 
cases there exists an implicit equation for the radius. A subsection deals with 
the simulation quality as well as along many directions. 



3.1 Spherical Symmetry - Conservation of Momentum 



The thin layer approximation assumes that all the swept-up gas accumu- 
lates in an infinitely thin shell just behind the shock front. The conservation 
of radial momentum at distance R,P(R), requires that 



P(R) = Po(Ro) 



(2) 



where Rq is the initial radius and Pq(Rq) the radial momentum at Ro evaluated 
when t = to. The conservation of the momentum gives 

^nR 3 pR = P (Ro) , (3) 



where R is the velocity of the advancing shock, see lDyson. J. E. and Williams. D. A 



(119971 ) ; iPadmanabhanl (120011 ) and p is the density of the ambient medium. The 



law of motion is 



R = Ro(l+^(t-t )j . (4) 
and the velocity 



R = R U + 4^(t-t )\ * . (5) 



We can derive Rq from equation (4) and insert it in equation (5) 



The astrophysical units are: t± and £o,4 which are t and to expressed in 10 4 yr 
units, R pc and Ro tPC which are R and R expressed in pc, Rkms and Ro 7 kms which 
are R and Rq expressed in km/s. Therefore the previous formula becomes 



R kms - 24.49 - 1 R tc Kpc ( 1 + Kc Kpc \ 4 (7) 



An interesting quantity can be the swept mass at a distance R from the origin 
4 

M = -nR 3 p . (8) 



The density can be 



p = lA ni m H -^- = 1.4 x 1.6610- 2 V^- , (9) 

cm 6 cm 6 



where ri\ is the number density expressed in units of particle cm 3 and the 
factor 1.4 takes into account elements heavier than hydrogen. Introducing the 
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solar mass, M , the swept-up mass is 
M — 0.143 R 1 3 n 1 M e , 



(10) 



where R\ is the radius expressed in pc units. 



3.2 Exponentially varying medium 



The number density of an exponentially varying medium is described by 



where z is the distance from the equatorial plane, no is the number of particles 
at R = R and h is the scale height. 

A three-dimensional (3D) expansion will be characterized by the following 
properties 

• Dependence from the instantaneous radius of the shell on the latitude angle 
9 which has a range [-90° <-> +90°]. 

• Independence of the instantaneous radius of the shell from 0, the azimuthal 
angle in the x-y plane, which has a range [0° -H- 360°] . 

The mass, M, swept-up along a solid angle A Q, between and R is 



n(z) = n exp 



z 



(11) 



h 



M(R) = ——lAm H n I m (R) + -irR 3 Q n m H 1.4 



(12) 



where is the mass of hydrogen and 




R 



(13) 



where R is the initial radius. The integral is 



h (2 h 2 + 2 R h sin (0) + R 2 (sin (9)) 2 ) e 



Rq sin(fl) 
h 



Im{R) 




(sin(0)) 3 



(14) 
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Fig. 2. Decimal logarithm of swept-up mass when R = 0.05 pc, n\ = 1 and A = 1. 
Physical parameters as in Table 1. 



In this case, the swept mass is a function of the latitude angle 9 and can be 
plotted assuming A Q, — 1, see Figure 2. 

This numerical evaluation gives a simple explanation for the asymmetry 
of the nebula around ^-Carinae: a smaller mass being swept means a greater 
velocity of the advancing radius of the nebula. 

Conservation of momentum gives 



M(R)R = M(R )Rq 



(15) 



In this differential equation of the first order in R, the variables can be sepa- 
rated and an integration term-by-term gives 



R 



J M{r)dr = M{Ro)Ro x (t - t ) , 
Ro 



(16) 



where t is the time and t o is the time at F tp. The resulting non-linear equation 
J-jvx, is equation (17) in IZaninettil ( 120091 ). The radius R pc as function of the 
time, is found by a numerical method, in the case of 77-Carinae the age is 
158 yr and therefore i 4 — t ,4 — to/(10 4 ?/r)=0.0158. 
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3.3 Power law medium 



A possible form for a power law profile of the medium surrounding the 
Homunculus nebula is 

n(z)=n (^y a , (17) 



where z = R x sin(#) is the distance from the equatorial plane, R is the 
instantaneous radius of expansion, no is the number of particles at R = Rq 
and a is a coefficient > 0. 

The swept-up mass, M, along a solid angle A f2 between and R is 
M(R) = ^lAm H n I m (R) + ^nR 3 n m H 1.4 , (18) 

where 



where Rq is the initial radius. Integrating gives: 



^3 f Rsin(0) 



I m (R)= \ no ' . (20) 

The resulting non-linear equation J^l of motion expressed in astrophysical 
units can obtained by eqns. (15) and (16) and is 



J~nl — ~ Ro,pc 4 & 2 + Rpc Ro,pc"oi 2 + R pc Ro, pc 3 (sin (#)) a a 



+7.0 R ,pc 4 a - Ro, pc 4 (sin (6)) a a-7.0 R pc R 



pc 



a 



+3 R 0iPC A (sin (9)y a - 12 R 0iPC 4 - 4.0 R pc R J (sin (0))- Q 



+ R pc - a+4 -°R ,pc a (sin (#))-" + 12 R pc R 



pc 



—0.122 Ro,p C 3 Ro,kms (t4 ~ h, 4) 
—0.01 Ro,pc 3 Ro,kms (*4 — ^,4) 

+0.0714 R ,pc 3 Ro,kms (U - t , 4 )a = , (21) 

where £4 and £q,4 are £ an d to expressed in 10 4 yr units, R pc and -Ro.pe are R 
and -R expressed in pc, Rkms and Ro,kms are i? and i? expressed in km/ s and 
is expressed in radians. 
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It is not possible to find R pc analytically and a numerical method must 
be implemented. In our case, in orde r to find the root of J-'nl, the FOR- 



TRAN SUBROUTINE ZRIDDR from lPress et all (I19921 ) has been used. The 



unknown parameters, Ro )PC and Ro t kms, are found from different runs of the 
code, t± — to 4 is an input parameter. 



3.4 Quality of the simulations 



As in IZaninettil (120091 ). two parameters are introduced to assess the quality 



of the simulations. The first, e, compares observed and simulated quantities: 
e= (1- l(^° bs ~^ um )l ). 1 oo ; (22) 

Rpc,obs 



where R pCt0 bs is observed radius, in parsec and -R pCi num is the radius from our 
simulation in parsec. 

The second defines an observational reliability, e Q b s , over the whole range 
of the latitude angle 9, 



lnn /-, J2j\Rpc,obs Rpc,num\j 

Cobs - J-UU^l — — 

2_,j -ftpc,obs j 



(23) 



where the index j varies from 1 to the number of availabl e obser vations. Those 



of the Homunculus are represented by Table 1 in ISmithl (120061 ) . 



3.5 Simulations quality for exponentially varying medium 



A typical set of parameters which allows the Homunculus nebula around r\- 
Carinae to be simulated in the presence of a medium whose density decreases 
exponentially is reported in Table 1. Table 2 presents numbers concerning the 
quality of fit. 

The bipolar character of the Homunculus is shown in Figure 3. In order 
to better visualize the two lobes, Figure 4 and Figure 5 show the radius and 
velocity as a function of the angular position 9. The ori entation of the observe r 



is characterized by the three Euler angles ($, 0, \l/), see lGoldstein et al.l ( 120021 ); 
different Euler angles produce different observed shapes. 

The velocity field is shown in Figure 6. 
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Fig. 3. Simulations lead to this picture of the Homunculus for an exponentially 
varying medium. The orientation of the figure is characterized by the Euler angles 
, which are $=130°, 0=40° and \l/=-140 o . Physical parameters as in Table 1. 




6 (degree) 



Fig. 4. Radius as a function of latitude for an exponentially varying medium (dotted 
line) and astronomical data with error bar. Physical parameters as in Table 1. 
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Table 1 

Parameter values used to simulate the observations of the Homunculus nebula for 
a medium varying exponentially (first 4 values) or a power law (2nd set of 4 values) 



Initial expansion velocity, i?o,i [km s x ] 8000 

Age (t 4 - to A ) [10 4 yr] 0.0158 

Scaling h [pc] 0.0018 

Initial radius Rq [pc] 0.001 

Initial expansion velocity, i?o,i [km s _1 ] 40, 000 

Age (t 4 - io )4 ) [10 4 yr] 0.0158 

Initial radius Rq [pc] 0.0002 

Power law coefficient a 2.4 



Table 2 

Agreement between observations and simulations for the Homunculus nebula, for 
an exponentially varying medium. 



radius velocity 
e(%) — polar direction 97 99 

e(%) — equatorial direction 87 19 
e obs (%) 85 75 



The accuracy with which our code reproduces the spatial shape and the 
velocity field over 18 directions of the Homunculus nebula as given by for- 
mula (23) is reported in Table 2. From a careful analysis of Table 2 we can 
conclude that the spatial shape over 18 directions is well modeled by an ex- 
ponential medium , e Q b s = 85%. The overall efficiency of the field is smaller 
tobs = 75%. We can therefore conclude that formula (23) which gives the ef- 
ficiency over all the range of polar angles represents a better way to describe 
the results in respect to the efficiency in a single direction as given by formula 
(22). 

3. 6 Simulations quality for power law medium 

For assumed parameters see Table 1, Table 3 reports the accuracy of radius 
and velocity in two directions. 
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Fig. 5. Velocity as a function of latitude for an exponentially varying medium 
(dotted line) and astronomical data with error bar Physical parameters as in Table 1 . 

Table 3 

Agreement between model for a power law medium and observations. 



radius velocity 
e(%) — polar direction 78 76 

e(%) — equatorial direction 6 2 
e ob s{%) 79 73 



4 Image 



An image of an astrophysical object is composed in a simulation by com- 
bining the intensities that characterize different points. For an optically thin 
medium the transfer equation provides the emissivity to be multiplied with 
the distance in the line of sig ht. The transfer equation has been analyzed in 



Section 5.1 of lZaninettil (120091 ) . The Homunculu s nebula was observed through 



emission-line spectra such as H 2 and [Fen], see ISmithl ( 120061 ). We now outline 
a possible source of radiation. The volume emission coefficient of the transition 
321 is 

n 2 A 2l hv 2l 

to = —tor- . (24) 



where level 1 is the lower level, level 2 is the upper level, n 2 is the gas number 
density, n 2 A 2 i is the rate of emission of photons from a unit volume, A 21 is 
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642.8433 km/sec 




S 18.3652 km/sec 

Fig. 6. Map of the expansion velocity for an exponentially varying medium. Physical 
parameters as in Table 1. 

the Einstein coefficient for the transition , h is the Planc k constant and V21 
is the frequency under consideration, see iHartiganl (120081 ) . In the case of an 
optically thin medium, the intensity of the emission 7 2 i is the integral along 
the line of sight 



'21 



hxdl 



(25) 



In the case of a constant gas number density 
hi oc I , 



(26) 



where I is the appropriate length, which in astrophysical diffuse objects de- 
pends on the orientation of the observer. 
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The numerical algorithm which allows us to build the image is now outlined. 



An empty (value=0) memory grid k) which contains NDIM 3 pixels 

is considered 

We first generate an internal 3D surface by rotating the ideal image 180° 
around the polar direction and a second external surface at a fixed distance 
AR from the first surface. As an example, we fixe d AR = 0.03 - R mrtr , where 
Rmax is the maximum radius of expansion, see ISmith et al.l (120061 ) . The 
points on the memory grid which lie between the internal and external 
surfaces are memorized on k) (value=l). 

Each point of k) has spatial coordinates x,y, z which can be repre- 

sented by the following 1x3 matrix, A, 



A 



x 



(27) 



The orientation of the object is characterized by th e Euler angles ($, 6, ^ 
and therefore by a total 3x3 rotation matrix, E, see iGoldstein et al.l ( 120021 ) 
The matrix point is represented by the following 1x3 matrix, B, 



B = E ■ A 



(2* 



• The intensity map is obtained by summing the points of the rotated images 
along a particular direction. 

An ideal image of the Homunculus nebula having the polar axis aligned 
with the z-direction which means polar axis along the z-direction, is shown in 
Figure 7 and this should be compared with the H2 emission structure reported 
in Figure 4 of ISmith! (120061 ). A model for a realistically rotate d Homunculus is 



shown in Figure 8. T his should be compared with Figure 1 in 
(l2000h or Figure 1 in lSmithl (120061 ). 



Smith fc Gehrz 



The rotated image exhibits a double ring and an intensity enhan cement 



i n the central region which c haracte r izes t h e little Homunculus, s ee ISmith 



(120021 ); llshibashi et all (120031 ): ISmithl (120051 ): iGonzalez et al.l (120061 ). Figure 9 
and Figure 10 show two cuts through the Homunculus nebula without and 
with rotation. The intensity enhancement is due to a projection effect and is 
an alternative for the theory that associates the little Homu nculus with an 



erupti o n occu r ring s ome time after the Great Eruption, see llshibashi et al 



(120031 ); ISmithl (120051 ) . We briefly recall that a central enhancement is visible 
in one of the various morphologies characterizing planetary nebulae. This can 
be compared with the model BL\ — F in Figure 3 of the Atlas of synthetic 
line profiles by iMorisset fc Stasinskal (120081 ) . 
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Fig. 7. Map of the theoretical intensity of the Homunculus nebula in the presence 
of an exponentially varying medium. Physical parameters as in Table 1. The three 
Euler angles characterizing the orientation are $=180°, 0=90° and ^=0°. This 
combination of Euler angles corresponds to the rotated image with the polar axis 
along the z-axis. 

Such cuts are common when analyzing planetary nebulae. As an example 



Figure 4 in iJacoby et al.l (120011 ) reports a nearly symmetrical profile of the 



intensity in the [OIII] image of A39, a nearly spherical planetary nebula. An- 
other example is the east- west cut in Hf3 f or the elliptical Ring ne b ula, c rossing 



the center of the nebula, see Figure 1 in iGarnett &: Dinersteinl ( 120011 ). Such 
intensity cuts are not yet available for //-Carinae and therefore can represent 
a new target for the observers. 



5 Conclusions 



Law of motion We have analyzed the law of motion determined by the 
conservation of radial momentum under two circumstances: exponential or 
power law variation of the density of the interstellar m edium. The two non- 
linear equations of motion are given by formulae (17) in IZaninettil f 20091 ) and 
21). A comparison of the observed and simulated profiles (Tables 2 and 3), 
shows that the profiles are better described by an exponential variation than 
by a power law. Here we have assumed that all the swept-up mass resides 
in a thin shell beyond the advancing surface. In the case of clumpy material 
this approximation can be refined by introducing the porosity p. Then, the 
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pc from Nebula center 

Fig. 8. Model map of the Homunculus nebula rotated in accordance with the ob- 
servations, for an exponentially varying medium. Physical parameters as in Table 1. 
The three Euler angles characterizing the orientation of the observer are $=130°, 
0=40° and ^=-140°. This combination of Euler angles corresponds to the observed 
image. 

I I I I I I I I I I I I I I I I I I I I I I I I I I I — i— I— I— I— I— I— I— I— I— I— I— r 




R (pc) 

Fig. 9. Two cuts of the model intensity across the center of the Homunculus nebula 
for an exponentially varying medium: equatorial cut (full line) and polar cut (dotted 
line). Parameters as in Figure 7. 
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Fig. 10. Two cuts of the model intensity across the center of the realistically rotated 
Homunculus nebula for an exponentially varying medium: equatorial cut (full line) 
and polar cut (dotted line). Parameters as in Figure 8. 

accumulated mass is M x l p rather than M. 



We briefly recall that the general principles of mass additio n to astrophysi- 
cal flo ws vi a hydrodynam ic mixing have been investigated by lHartquist et al 
(119861 ) and iPittardl (120071 ). In order to fix the porosity parameter p the equa- 
tion of motion along the polar direction should be provided in the form 



R{t) = r obs t 



Oiobs 



(29) 



where the two parameters r b s and a a b s are found from an analysis of the 
observational data. As an example the data over a 10 year period of SN 1993 J 
can be approximated by a power law dependence of the type R oc t 0,82 , see 
Marcaide etldl (120091 ). 



Images Assuming an optically thin medium, it is possible to make a model 
image of the Homunculus nebula once two hypotheses are made: 

(1) The thickness of the emitting layer, AR, is the same everywhere AR = 
0.03R max , where R max is the maximum radius of expansion. 

(2) The density of the emitting layer is constant everywhere 

A 2D image of the Homunculus nebula is shown in Figure 7 and a non-rotated 
image in Figure 8. Our model provides an explanation for the emission maps 
in e.g. X-rays ( 0.2 -1. 5 keV) by CHANDRA. An inspection of Figure 1 in 
Corcoran et al.l (J2004aj) shows that the next generation of X-ray observatories 
may resolve the double ring as observed in the near infrared, see Figure 4 of 
Smith! (120061 ) and our simulation in Figure 8. The validity of our assumptions 
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can be checked by looking at cuts in intensity over the various bands. We 
recall that cuts in intensity allow to fix the thicknes s of the swept-up la yer of 
a planetary nebula such as A39 where^ ~ -^p, see IJacoby et al.l (120011 ) . 
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